Influence of the glass transition on the liquid-gas spinodal decomposition 
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We use large-scale molecular dynamics simulations to study the kinetics of the liquid-gas phase 
separation if the temperature is lowered across the glass transition of the dense phase. We observe a 
gradual change from phase separated systems at high temperatures to nonequilibrium, gel-like struc- 
tures that evolve very slowly at low temperatures. The microscopic mechanisms responsible for the 
coarsening strongly depend on temperature, and change from diffusive motion at high temperature 
to a strongly intermittent, heterogeneous and thermally activated dynamics at low temperature, 
leading to logarithmically slow growth of the typical domain size. 
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The kinetics of phase separation is a fascinating sub- 
ject with a long history in condensed matter physics 
Although this nonequilibrium phenomenon is complex, 
universal scaling behaviour have been indentified and can 
be studied using statistical mechanics (H. Phase separa- 
tion kinetics also has important experimental and tech- 
nological consequences [3j. While this phenomenon is 
well-understood for simple materials, more complex be- 
haviour can be expected if the different phases have spe- 
cific properties, such as viscolelastic polymer solutions or 
colloidal suspensions, which interfere and possibly mod- 
ify the coarsening process in a non-trivial manner. 

Here we study numerically the kinetics of the liquid-gas 
spinodal decomposition over a broad temperature range 
encompassing the glass transition of the dense fluid. Not 
much is known about the resulting gas-glass spinodal 
decomposition-a specific instance of a viscoelastic phase 
separation although this interplay is invoked in a 
number of contexts @, H. Further interesting cases are 
colloids Q and proteinsj^ with short-ranged attraction 
that form a gel, instead of phase separating. A plau- 
sible interpretation is that the dense phase is actually a 
glass, and this 'arrests' the phase separation 0,|f|. It has 
recently been suggested that similar structures could be 
obtained in Cgo molecules 0], and that the range of the 
potential is important Q. Bicontinous disordered struc- 
tures reminiscent of the ones obtained in spinodal decom- 
positions may also be found in colourful bird feathers, 
and were recently interpreted as incompletely phase sep- 
arated polymeric glasses Q . Although it can be expected 
that the glass transition slows down the spinodal decom- 
position, the general mechanisms at play have never been 
studied in any detail on the microscopic scale. 

Although coarsening processes are most efficiently 
studied theoretically using coarse-grained models, such 
as model H for the liquid-gas spinodal decomposi- 
tion [l[, in practice it is difficult to faithfully incorpo- 
rate the complex — typically highly nonlinear and history 
dependent — physical properties of glasses in such descrip- 
tions [To}. Thus, we work in the opposite direction, 
and start from a microscopically realistic description of 



the homogeneous glass and study its behaviour during 
phase separation. Indeed, numerous successful atomic- 
scale simulations of the liquid-gas spinodal decomposi- 
tion exist 11 -Til- In Ref. [HI], a Lennard- Jones system 
was quenched to low temperature in the coexistence re- 
gion, and the resulting gas-crystal phase separation was 
studied, but no arrest was reported. Simulations of re- 
alistic colloidal interactions were reported 
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the quenches have been performed at very low temper- 
atures where particle aggregation are nearly irreversible 
and thermal fluctuations play little role. 

We use molecular dynamics simulations to study a 
80:20 binary Lennard- Jones mixture with interaction pa- 
rameters chosen to yield excellent glass-forming abil- 
ity [lj|. Specific attention was paid here to system sizes. 
While numerical studies of the glass transition in the ho- 
mogeneous liquid typically require simulating about 10 3 
particles, we found that up to 10 6 particles were needed 
to obtain results devoid of finite size effects during the 
phase separation. We obtained most of our quantitative 
results using 3 • 10 5 particles. We study the kinetics of 
the phase separation following an instantaneous quench 
at constant density from the high temperature fluid phase 
to the coexistence region, which we determined for the 
present system. In the following we use Lennard- Jones 
units corresponding to the majority component, express- 
ing length in units of particle diameter, a, and time in 
units of r = y/ ma 2 /e, where m is the particle mass, and 
e the energy scale in the Lennard- Jones interaction. 

In Fig. Q] we present a density, p, and temperature, T, 
phase diagram showing the location of the homogeneous 
gas, liquid, and glass phases of the system, as well as the 
coexistence region. For the glass 'transition', we show 
the location of the temperature T mct obtained numeri- 
cally [i^l from fitting the temperature evolution of the 
relaxation time of the system using predictions from the 
mode-coupling theory. Usually, it is hard to maintain 
thermal equilibrium below T mct , which can be taken as 
a proxy for the 'computer' glass temperature. Figure [1] 
demonstrates that the present model exhibits the right 
ingredients, as in a quench below T m 0.3 in the coexis- 
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FIG. 1: Phase diagram of the binary Lennard- Jones mixture 
showing the location of the homogeneous gas, liquid, and glass 
phases, and the gas-liquid and gas-glass coexistence regions. 
The glass line corresponds to the mode-coupling temperature 
(Ref. [U), while the spinodal is taken from Ref. 0. The 
vertical dashed lines separate disconnected droplets from bi- 
continuous structures. 



tence region, the dense phase is a glass — not a liquid. 

On the liquid side, the coexistence line hits the glass 
line near (p « 1.15, T w 0.3). The glass line can be 
followed between the binodal and the spinodal, where 
the homogeneous liquid is metastable. Using methods 
described below, we found that the coexistence line is 
only weakly affected by the glass line, as found in some 
experiment Wc do not find a reentrant coexistence 
line slaved to the glass line, as reported in Ref. 

In Fig. [5] we present typical configurations obtained in 
various parts of the phase diagram. Figure Ot shows a 
homogeneous configuration obtained in the liquid phase, 
while Fig. [5b shows the final configuration obtained after 
a quench to (T — 0.5, p — 0.4), where the system is fully 
demixed at long times. It is in this temperature regime 
that the liquid-gas spinodal decomposition is typically 
studied [ll-[ll|. For a quench to densities below p m 0.2 



we find that the coarsening proceeds via growth of dis- 
connected droplets, as illustrated in Fig. [2J;. Similarly, 
we find that for densities above p « 1.0, the coarsen- 
ing proceeds via the growth of disconnected gas bubbles 
immersed in a dense fluid. More interesting is the low 
temperature regime in the density range p G [0.2, 0.8] 
where bicontinuous, percolating, gel-like structures are 
obtained at long times, see Fig. This means that 
phase separation is incomplete when T is low, at least in 
the timescale of the simulation. These configurations are 
strikingly reminiscent of the confocal microscopy images 
in colloidal gels 0, |f| . 

To study quantitatively the phase separation kinetics, 
we must characterize the structures shown in Fig. [2ji and 
follow their time dependence for a quench at a given state 
point. Phase separating systems are usually character- 




FIG. 2: Representative structures of (a) the liquid phase 
at (T = 1.75, p — 0.6), (b) a fully demixed system at 
(T = 0.5, p = 0.4), (c) disconnected growing droplets at 
(T = 0.1, p = 0.1, t — 10 4 ) (d) a bicontinous gel-like struc- 
ture obtained at long time at (T = 0.1, p = 0.2, t = 10 4 ) 
where phase separation is incomplete. 



ized in numerical work using the static structure factor 
of density fluctuations, or equivalently, the pair corre- 
lation function. However, we found that a numerically 
more efficient characterization of the bicontinuous struc- 
tures is provided by the chord length distribution of the 
low-density phase, which contains detailed information 
about the domain size distribution. This is a standard 
characterization tool for porous media (22. 23 1 . 

To measure this quantity, we must first numerically de- 
termine the location of the interfaces separating the two 
phases. This requires a non-trivial coarse-graining for 
a particle-based simulation in which complex topologies 
are present. Thus, we discretize space in cubes of small 
linear size d = 0.5a, and define on this cubic lattice a 
coarse-grained density field, p(r), using a weighted aver- 
age over cells around r. The distribution of the coarse- 
grained density, p(p), is typically bimodal, as expected 
for a two-phase system. In particular, p(p) exhibits at 
large p a maximum, which provides an accurate measure 
of the average density of the dense phase. We have used 
this maximum to determine the coexistence line at low 
temperatures in the phase diagram of Fig. [1] Moreover, 
a careful analysis of p{p), backed by direct visualizations, 
shows that a density threshold can be estimated to accu- 
rately locate gas and fluid domains, and the position of 
the interfaces separating them. 

The time evolution of the so obtained interfaces is 
shown in Fig. [3] after a quench at p = 0.4 in the coex- 
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FIG. 3: Time series for a quench at p = 0.4 showing interfaces, 
the yellow/green side being inside the gas/fluid phase. While 
the liquid-gas spinodal decomposition at T — 0.5 proceeds 
until complete phase separation, dynamics is strongly slowed 
down in the gas- glass case at T = 0.1. 

istence region. At T = 0.5, the spinodal decomposition 
proceeds rapidly until complete phase separation of the 
system. By contrast, for T = 0.1 the evolution is similar 
to the one at T = 0.5 at short times, t < 10 2 , but further 
evolution is strongly suppressed (but not fully arrested) 
over the remaining two decades of the simulations, up to 
t = 10 4 . 

We are now in a position to measure the chord length 




t 

FIG. 4: Evolution of the typical domain size at p = 0.6 and 
various T, in a log-log plot. Below T ~ 0.3, the growth is 
much slower than the diffusive law, L(t) ~ V s (dashed line), 
but not fully arrested. 

distribution and its evolution. Chords arc defined by two 
consecutives intersections of a straight line with the in- 
terfaces of the two-phase material. In practice, we mea- 
sure chords along the three axis of the lattice used to 
coarse-grain the density, and measure the length £ of 
the segments belonging to the gas phase. Thus, we ob- 
tain for each configuration a distribution P(£,t). We 
find that P(£, t) has a maximum, and beyond the max- 
imum is well-described by an exponential tail, as com- 
monly found in porous media [H, H3|- We define the 
average domain size as the first moment of the distribu- 
tion: L{t) = f o °° £P(£,t)d£. Wc find that L(t) is a robust 
and efficient measure of the typical domain size, which 
could be easily implemented in confocal microscopy ex- 
periments ys| . 

In Fig.HJwe report the time evolution of L(t) at density 
p = 0.6 and various temperatures from T = 0.6 down to 
very low temperatures, T = 0.015. When temperature 
is large, we find that the growth law is well described 
by a power-law, L(t) ~ t n , with an apparent exponent 
n ~ | (T1I - IT3I . We interpret this behaviour as the indi- 
cation that our data are taken in the crossover between 
the diffusive (n = -|) and hydrodynamic (n = 1) regimes, 
which arc notoriously difficult to disentangle in molecular 
dynamics simulations [3, HH . 

The most important observation in Fig. |4] is that for 
T < 0.3 the growth law is always much slower than the 
diffusive t 1 ^ 3 law and strongly depends on temperature, 
as reported in a recent colloidal experiment Q . The data 
are actually curved in a log-log representation, suggest- 
ing that the growth is not algebraic, but presumably log- 
arithmic. We emphasize that these data are free of finite 
size effects. Indeed we found that smaller system sizes 
yield asymptotic lengthscales that are smaller, and spin- 
odal decomposition appears to be much more strongly 
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(a) T=0.5 



FIG. 5: Mobile particles (green spheres) and their displace- 
ments (red cylinders) at p = 0.4 for time intervals correspond- 
ing to a equivalent domain sizes, (a) T = 0.5, times t = 32 
and t = 34.5. (b) T = 0.1, times t = 1096 and t = 1905. 



slowed down at low T if N is chosen too small, a fact 
which is experimentally relevant for gelation in confine- 
25| . A closer look at the data at large times and 
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low T in Fig. [4] shows that if the growth is much slower 
than diffusive, it is not completely suppressed either, that 
is, spinodal decomposition is not kinetically arrested by 
the glass transition of the dense phase, as assumed in pre- 
vious work 0-0, [HI , although it does become difficult to 
detect at very low temperatures. 

Instead, we find that low temperatures, T < 0.3, 
strongly modify the microscopic mechanisms of the phase 
separation. While for the liquid-gas phase separation 
particles diffuse to relax the curvature of the interfaces, 
at low T surface tension seems to play very little role. 
Direct visualization shows that curved interfaces quickly 
solidify and hardly relax, while particle diffusion in the 
dense phase (an aging glass) is strongly suppressed. 

Particle displacements help to understand the micro- 
scopic mechanisms at play. In a given time interval, we 
select about 1 % of the particles that have moved far- 
thest. At high T we find displacements as in Fig. [5^ 
where mobile particles appear to be randomly scattered 
throughout the system, with uncorrelated displacements. 
The picture is very different at low T, see Fig. [5b, which 
reveals strong clustering of the fast moving particles, 
with very correlated displacements. At low temperature, 
particles in the dense glassy domains are basically ar- 
rested and can only move significantly if the whole do- 
main moves as a rigid body. Such large displacements 
occur in a temporally intermittent and spatially heteroge- 
neous manner, and follow a change in the local topology 
of the bicontinuous structure, for instance when a thin 
neck breaks. There are three such relaxation events in 
the example of Fig. [5b. We speculate that the main driv- 
ing force of the phase separation process is the mechan- 
ical stress stored in these non-equilibrium, disordered, 
porous structures, rather than surface tension. Due to 
thermal fluctuations, these constraints are released in a 
highly heterogeneous manner, yielding large scale rear- 



rangements with correlated particle displacements, which 
ultimately increase the average domain size. 

In future work, one should characterize in more detail 
the microscopic aging dynamics of the system and extend 
our studies to different kinds of models, as spinodal de- 
composition represents a promising tool to produce dis- 
ordered porous media, even in atomic systems. 
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